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Abstract. We reconstruct the time dependent quantum map corresponding to the 

relaxahon process of a two-spin system in liquid-state NMR at room temperature. 

By means of quantum tomography techniques that handle informational incomplete 

data, we show how to properly post-process and normalize the measurements data 

for the simulation of quanfum informahon processing, overcoming the unknown 

number of molecules prepared in a non-equilibrium magnetisation state (Nj) by an 

initial sequence of radiofrequency (RF) pulses. From the reconstructed quantum map, 

we infer both longitudinal (Ti) and transversal (T2) relaxation times, and introduce 

the /-coupling relaxation times (t/ T^), which are relevant for quantum information 

processing simulations. We show that the map associated to the relaxation process 

carmot be assumed approximated unital and trace-preserving for times greater than 

jJ 

l2- 
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1. Introduction 


During the last few years, we have seen proposals and experimental implementations 


of quantum information processors for various quantum systems I13I. While there 
is still no ultimate physical platform, nor an underlying information processing 
model, most of the known candidates must deal with adverse quantum effects. 
Nuclear Magnetic Resonance (NMR) techniques have been successfully used to 


implement quantum protocols, e.g., quantum teleportation and Shor's | [^ 
algorithm. NMR remains indispensable in quantum computing as it continues to 
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provide new insights, new methods and new techniques, as well as allowing for 
testing interesting quantum computing tasks, although there are still huge difficulties 
in initialization, error-correction and scalability |^ . 

In the NMR Ensemble Quantum Computing (EQC) introduced by Cory et al. Q, 
a key observation is that, compared to the ideal case where one hundred percent of the 
molecules would be polarized, the proportion of molecules effectively contributing to 
the magnetisation in a pseudo-pure state is about one part in a million, nonetheless 
one still has a macroscopic detectable signal. On the other hand, when performing 
an algorithm with EQC, one has to deal with two intrinsic limitations of the liquid 
state NMR, namely: (z) the natural relaxation process the spins are subjected to 
is non-unital, which means that the state simulated by a pseudo-pure state has a 
non-linear time evolution, what can compromise the quantum computation, (z’z) 
Related to this fact, there is the recovery of the magnetisation towards equilibrium, 
characterized by the spin-lattice relaxation time Ti, which diminishes the amplitude 
of the component corresponding to the simulated state along the time. This means 
that expectation values measured in different rounds of the experiment will not 
be comparable, xmless one takes into account the reduction of their magnitudes 
along the time. When the quantum computation lasts about the transverse relaxation 
time T2, the typical coherence time, these drawbacks are negligible, and one can 
ignore both the non-linear evolution and the degradation of the simulated state, 
and also some particular expectation value can be taken as a reference in order 
to normalize the other measurements in the computation. Notwithstanding, the 
spins natural relaxation is an ever present process relax) superimposed to the map 
of the particular computation one is performing. In the present work, our goal is 
to characterize the relaxation process in the same way one would reconstruct the 
map corresponding to a quantum computation, i.e., by means of quantum process 
tomography with the least possible assumptions. We will be able to quantify, along 
the time, how much the relaxation process departs from being unital and trace 
preserving, and we will also characterize how the relaxation process affects the 
normalization of expectation values. This way we will learn in which time window it 
is safe to assume the relaxation process is approximated unital or trace preserving, 
complementing other previous studies |j3 ^ Thus we hope to shed some light in 
how one could compensate for the two drawbacks aforementioned, and consequently 
increase the time window for quantum simulations in liquid state NMR. The present 
work also serves the purpose to test, in a real NMR experiment, and in a non trivial 
situation, the efficiency, robustness, flexibility and simplicity of implementation of 
our quantum tomography techniques 

Robust techniques for quantum process tomography in NMR EQC have been 
developed in the case the process is assumed to be Markovian and trace preserving 
[|^ , or when some model is assumed for the Hamiltonian governing the process [loj. 
In both cases, the process tomography was performed in a time window of about T2, 
and therefore the process could be safely assumed to be unital. We, on the contrary. 
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will assume the worst case scenario, meaning the process may be non-unital, non- 
Markovian, and not trace preserving. Besides that, we want to follow the quantum 
process for very long times, much more than the spin-lattice relaxation time Ti, what 
raises the problem of how to compare measurements taken in different times, which 
is what we call the normalization problem. A pedagogical discussion of why the state 
simulated by a pseudo-pure state evolves non-lienarly under the relaxation process, 
and an ad hoc phenomenological method to circumvent the normalization problem 
were presented by Gavini-Viana et al. |j^. The basic idea is that a pseudo-pure state, 
according to Cory et al. has the form 

Ppp = [{ 1 -a)l 2 n+ 2 a\xp){xp\]/[{l-a) 2 ^+ 2 a], (i) 

where n is the number of spin-^ particles and 12” is the 2” x 2” identity matrix]^ In 
a quantum computing simulation, we are concerned with unitary evolutions, which 
of course affect only \ip){xp\, the simulated state. However, if 0 is a non-unital map, 
‘h(l2«) 7^ I2"/ rewriting (^{ppp) as a pseudo-pure state is equivalent to a non-linear 
evolution of \ip){xp\. Note also that the new pseudo-pure state ^{ppp), differently 
from the case of unital evolutions, has a new "a". Therefore, the normalization problem 
is equivalent to account for the time evolution of "oi” under non-unital processes, 
such that we can always track the evolution of the resulting pseudo-pure state. 

To avoid confusion, it is important to remember that a quantum process is a 
complete positive map, that can be expressed in the usual Kraus sum representation, 
namely: 

<t>{p} = Y^KipKl (2) 

i 

where the K, are the Kraus operators. A map is unital if and only if 

= I 2 ., (3) 

i 

a condition that does not imply trace preservation. To be trace preserving, a map 
must satisfiy 

= I2.. (4) 


The state simulated by a pseudo-pure state evolves linearly under unital and non¬ 
trace preserving maps, and therefore these maps pose no particular challenge to 
usual process tomography techniques in liquid state NMR EQC. In particular, the 
quantum process tomography method developed by Maciel and Vianna |i^, and 


experimentally tested in Quantum Optics by Vianna et al. I30I, could be employed 
without any modification to handle these maps, as well as non-unital ones. A caveat 
about our basic assumption of the complete positiveness (CP) of the map is that the 
tomographed map could eventually be non-CP due to noisy measurements 1 31 1 , or in 


the case of correlations between the initial prepared state and the environment [261. 


f One should remember that the pseudo-pure state is not an actual quantum state of the NMR system, 
but rather the result of spatial or temporal labelling II22). 
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The methods we will employ here | [i^ i8 exploit the natural constraint 
of positive semidefiniteness arising in quantum tomography, in tandem with 
semidefinite programs (SDP) 0 , which are convex linear optimization problems with 
linear matrix constraints. A great merit of an SDP is that to find its global minimum 
there are many available implemented numerical methods which converge efficiently 
in polynomial time Ijq ^ 2^. Our tomographic methods [19 18 ^ are very easy 
to implement and are also very flexible. They consist of the minimization of some 
figure of merit, expressed as a linear objective function, subject to the constraint of 
the positive semidefiniteness of the density matrix, in the case of state tomography, or 
the positive semidefiniteness of the dynamical or Choi matrix, in the case of process 
tomography, with the restriction that the resulting state or process be compatible 
with the measurements. The choice of the objective linear function to be minimized 
makes our method equivalent to other approaches. For example, by choosing the 
nuclear (trace) norm of the positive semidefinite matrix of interest as the objective 
function, we have a matrix completion approach equivalent to compressed sensing 
IN. Note however, that nuclear norm minimization is also a good heuristic for 
quantum tomography in general, including the case of informationally complete 
measurements. Detailed comparison of our methods with maximum likelihood and 
maximum entropy was reported in Q. Our methods were tested for state and process 
tomography in Quantum Optics experiments ^6 25 30 1 , and in liquid state NMR it 
was tested for state tomography 1241 . 

The paper is organized as follows. In the next section, we describe the liquid 
state NMR setup we use to simulate two qubits, and the basic theory modelling our 
system. Section 3 is dedicated to a discussion of the normalization problem. Section 
4 contains our tomographic methodology. In Sec. 5, we present and discuss the 
results. We conclude in Sec. 6. 


2. The NMR system 

The NMR experiment was performed on a liquid-state enriched carbon-13 Chloroform 
sample, at room temperature, in a Varian 500 MHz shielded spectrometer. This 
sample exhibits two qubits encoded in the and spin-| nuclei. The 
heteronuclear interaction of the two spins is characterized by the magnitude of 
the /-coupling being much smaller than the resonant frequencies of the two nuclei, 
500 MHz and 125 MHz, for ^H and respectively. This system is useful for the 
simulation of quantum information processing in a Hilbert space of two qubits. 

The quantum model adopted in room temperature liquid-state NMR experiments 
is of an ensemble of N ( 0 ( 10 ^^)) non-interacting independent molecules. Therefore 
the Hamiltonian can be written as 

=H® + 1 O HO ^ _ 

... + (g) H o 1 + o H 
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N N 

^^{N} = E(8)(^y^^+(l-%)!)' (5) 

;=1 j=l 

where 1 is the 2 x 2 identity matrix, 3 ij is the Kronecker's delta, and H is an effective 
Hamiltonian for a single molecule. The wave function and the density matrix 
corresponding to the N-particle Hamiltonian are assumed to be separable, namely: 

N N 

'f{N} = 0 ^i and p^^^ = 0 pi, (6) 

i=l i=l 

where xpi and pi are the wave function and density matrix for a single molecule, 
respectively, where we shall work only in the spin degrees of freedom. 

In an NMR experiment, one does not have access to the particular state of a 
distinct molecule in the preparation procedure, all one has is the average state of a 
representative molecule of the ensemble. Therefore, we can write: 

N 

P{N} = ^pi = ( 7 ) 

i=l 


where p stands for this representative state, which we shall discuss further later. 

Our observables are related to the resultant macroscopic magnetisation induced 
by an external static magnetic field, and we assume they are operators of the form: 

= A O + 1 (g) A O + ... 

. . . + (g) A (g) 1 + (g) A, 

N N 

^{N} = E 0 + (1 - • (8) 

1=1 7=1 

Therefore, we have expectation values of the kind: 

An,) = Ef 

= N tr (|0A). 


(g){SijA, + {l-Sp)l) 


This rationale, assuming Eqs. (|5|j^ as premisses, formalises a model in which 
one has equally prepared, but independent, black boxes and the same measurement 
"buttons" are pressed in all of them when the measurement occurs. Thus, we have a 
reducible to a single box model and we can generalise and specialise our arguments 
back and forth, i.e., what happens in a single box, will happen in the other ones 
and vice versa. We shall use this freely. This model fits the common knowledge that 
measurements in NMR have contributions of all the molecules in the ensemble 1 22 ^ 
[2^ , justifying the detected macroscopic signals. 

The observables of interest in the NMR experiment are generally the components 
of the macroscopic nuclear magnetisation, which are proportional to the ensemble 
average values of the nuclear spin operators. The matrix representation of these 
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operators in the \S,m) basis (with S = 1/2) reads Sx = (Tx/ 2, Sy = o'y/2, Sz = Pz/2, 
with the usual Pauli matrices given by: 



1 

0 

1 


1 

0 

1 

T 


1 

1- 

0 

II 

1 

0 

II 

i 

0 

II 

1 

0 

-1 


For example, one can calculate the magnetisation in the z-direction of an ensemble of 
nuclear spins using the following relation: 

(Mz) = 7z oc tr {pcTz), (lo) 

where 4 is the empirical value associated with the observation of the magnetisation 
in the z-direction and p is the density matrix of the system of a single specimen. 
Remember that the left-hand side of this equation represents the statistical average 
over the entire ensemble of N molecules, and not the expectation value for a given 
particular member of the ensemble. 

The relevant effective two-spin Hamiltonian in Eq. is given by 

H = —hconSz 01 — hcoct ® Sz + InhJSz 0 Sz, (ii) 


where the labels H and C stand for hydrogen and carbon, respectively, as the first and 
second particles. The resonant frequencies oou and ooq include the effects of isotropic 
chemical shift Ufso in acetone solution for each nucleus, i.e., oo = col {1 — O'iso)- The 
Larmor frequency is given by ooi = 7 hBo/ where 7„ is the gyromagnetic ratio of the 
nucleus, and Bq is the static magnetic field in z-direction. 

At thermal equilibrium, the density matrix obtained from the Hamiltonian in Eq. 
( [ii] > is given by the maximum entropy principle |i^ : 

g-H/ksT 
= tr (g-H/ZcgT) ' 

where fcg is the Boltzmann constant and T is the temperature. In the high temperature 
limit, one can expand this exponential as 


Peq 






- Ihcoc-n 




(13) 


AReT ^RbT IReT 

Note that the last term of Eq. ( [i3] > is usually neglected (see for example {2^ 24 1 ), for 
it is actually very small compared to the other terms. Notwithstanding, as we will 
show, this term contains crucial information about the time window for quantum 
computing simulations. 

Now, we are able to calculate the total magnetisation of the system in the z- 
direction, using the equilibrium state in Eq. ( [i3] >, the relation in Eq. ( [1^ and Curie's 
law, which leads to 



^h7h^^^o I ^c7c^^Bo ^ njTZ^Ljh^Bo 
4 ReT ^ ARbT ^ R^ 


(14) 


where nu and uq are the number of hydrogen and carbon nuclei per unit volume, 
respectively, which are equal and correspond to the number of chloroform molecules 
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per unit volume {uh = = n); nj stands for the number of responsive pairs of 

/-coupled nuclei per unit volume, and lj its "gyromagnetic ratio". These responsive 
pairs of spins correspond to a small fraction of chloroform molecules {Nj/N 1 ), 
for which the spins of and are in a chosen coherent state prepared by a 
certain sequence of radiofrequency pulses (RF). The first issue on normalizing the 
experimental data arises due to the fact that nor Nj, neither its density Uj, are 
straightforwardly revealed by the magnetisation measurements. We shall discuss this 
later. 

Let us focus on the experiment now. We know that interesting NMR quantum 
information processing experiments happen in non-equilibrium, where the responsive 
/-coupled pairs of spins play an important role. In our experiment, we prepare the 
system in a desired non-equilibrium magnetisation state at time f = 0 and allow 
it to relax. First, at t = 0 , a certain sequence of RF pulses prepares the system 
in some desired non-equilibrium state, where we have the maximum number of 
responsive /-coupled spin pairs. This results in a small macroscopic non-equilibrium 
magnetisation. The second part of the experiment happens in the time window 
which we can use the small fraction of responsive molecules in the non-equilibrium 
state to simulate the dynamics of quantum parameters {e.g., I29I). The system is 
allowed to relax, the molecules that responded to the preparation start to recover 
their equilibrium state, diminishing Nj, and finally (f —^ 00 ) the system reaches the 
equilibrium with vanishing Nj. 

The relaxation is such that, after some time has elapsed, the magnetisation 
returns almost completely to the initial z-direction, satisfying again the thermal 
equilibrium requirements. It is worth emphasising that the exact meaning of "some 
time" is largely dependent on the details of each particular nuclear spin system 
and its environment, ranging typically from microseconds to several hours. The 
relaxation consists of two different processes, occurring simultaneously but, in 
general, independently, namely: the transverse relaxation and the longitudinal relaxation 
1 ^ . We shall characterize, by the quantum process tomography, a net effect of these 
two phenomena, that is paramount to the problem of NMR simulation of quantum 
information processing. For reasons to become clear soon, we shall refer to this net 
phenomenon as the J-coupling relaxation. 

The longitudinal relaxation is the process that leads the longitudinal component 
of the nuclear spin magnetisation to recover its equilibrium value. The recovery of 
the Mz component of the magnetisation is related to transitions between the nuclear 
spin levels. The natural tendency is the system to give up its excess of energy by 
effecting transitions from the upper to the lower energy level. After some time, which 
is commonly named Ti, the Boltzmann distribution is almost reestablished. 

Simultaneously to the longitudinal relaxation, the transverse relaxation is 
the process that leads to the disappearance of the components of the nuclear 
magnetisation M that are perpendicular to the static field Bq. The origin of the 
transverse relaxation relies on the loss of coherence in the precession motion of the 
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spins (or dephasing of the spins), caused by the existence of spread in precession 
frequencies tor the collection of nuclear spins. This spread progressively results in a 
reduction of the resultant transverse components (e.g., Mx and My). After some time, 
which is commonly named T2, the spins distribute randomly in a precession cone 
around Bq and the transverse magnetisation tends to zero. 

The /-coupling relaxation, which we introduce here, can be considered as the net 
effect of both longitudinal and transverse relaxations. In NMR quantum information 
processing simulation experiments, correlated spins in non-equilibrium states are 
indispensable to simulate quantum algorithms. We name r| and the characteristic 
times in which the coupled magnetisations Mij iy i,] G {x,y,z}) decay, for the 
longitudinal, or Mzz component, and for the transverse components. Note that 
what we call coupled-magnetisation refers specifically to expectation values of two-spin 
operators ^tr erf while uncoupled-magnetisation is related to expectation 

values of one-spin operators ^tr (perf (g) l) or tr ^|0l ^ erf Our motivation for 

introducing and was the following, (i) These parameters can be extracted from 
the usual NMR data; (ii) the Mzz component behaves distinctly from all the other 
coupled magnetisation components; (Hi) and finally, these parameters characterize an 
important change in the behaviour of the relaxation process map, as we will show. 

In order to reconstruct the density matrix of the system along the time, we 
need to normalize the experimental data, such that we can relate the measured 
non-equilibrium magnetisations to expectation values. Here, we are not concerned 
with the simulated states, we deal only with the real non-equilibrium NMR states, 
although it is worth to emphasise that, in the preparation procedure, we also assume 
ergodicity, as we collect the whole set of measurements using temporal labeling 

0 - 

3. The normalization problem 

The normalization issue is a daunting problem in NMR experiments. Firstly, by 
the data acquisition itself, since it involves taking the intensity of the peaks from a 
Fourier transform spectra of an oscillating signal. Secondly, because we are forced 
to assume the representative state p in Eq. Q. It would be preferable if we could 
address each molecule individually like 

P{N}it) = Pl{t)®p2it)...®pNit). 

If it were the case, we would be able to coimt how many spins are still coherent and 
responding accordingly, therefore it would be possible to know exactly the number 
of responding molecules as in Eq. (|^ {e-g-, if we were interested in knowing Nj). 
However such a case is not possible and we shall stick to the weaker premise of the 
representative state, losing track of this counting. Thirdly, as consequence of the 
traceless nature of the operators corresponding to the magnetisation components. 


122 24 
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we do not possess any direct experimentally measured data related to the trace of 
the non-equilibrium density matrix. If we were measuring the eigenprojectors (Pj^) of 
these operators (A = Ylk=i ^kPk)/ this problem would not appear, as the completeness 
relation = Yly. P^ would give us a trivial recipe for normalization, i.e., we would 
impose Yi:=i h {p^k) = 1 - However, what we have is (A) = tr (|0A), and for this 
measurement to be useful, we need some reference expectation value (Ayef). Of 
course, the reference expectation value cannot vary in time. For example, in the EQC 
MMR experiment reported in which employed the same experimental setup 


and molecule of the present work, the reference expectation value was O ) 
in a pseudo-pure state simulating the two-qubit state | tt) ■ It was imposed then 
that /^ref = 1/ other expectation values were divided by Nygf. 

As we shall see, the quantum tomography techniques we employ here ^ 
allow us to circumvent this problem, yielding proper inferences of the normalized 
non-equilibrium quantum states and quantum maps. 

We have three different kinds of signals corresponding to our measurements, 
which we have to deal with in order to normalize the experimental data. One may 
grasp this in Fig. We have, responding to RF pulses in the resonant frequency coh 
of the hydrogen, both coherently coupled and uncoherently coupled hydrogen nuclei; 
the same happens to carbon, for RF pulses in the resonant frequency coq- As the 
gyromagnetic ratio of the hydrogen is almost 4 times larger than that of the carbon, 
the hydrogen species is expected to have a larger magnetisation in equilibrium. The 
last kind of signal, which is very small compared to the first two, corresponds to that 
coherently coupled magnetisation of the Nj chloroform molecules. To remedy the lack 
of some relation tying the intensity of measured magnetisation components, like the 
completness relation we have mentioned in the previous paragraph, we normalize all 
measured signals I15I corresponding to non-equilibrium magnetisations backwards, 
such that the expected value of the magnetisation in equilibrium is equal to one. 
Therefore, the equilibrium magnetisation plays the role of our reference expectation 
value that does not vary in time, as discussed in the previous paragraph. This 
rescaling of the magnetisation signals is also convenient in order to perform our 
numerical optimizations with good precision. For the hydrogen and carbon signals, 
we have 


■H,: 


Pit) = 


1 hWH / ‘“'(t) 


fHeq 

J -7 


M4:kbT 

and for the coupled signals we have 


and 7 y'(f) = 


1 hcocljii) 


M 4 :kuT 


fCeq 


(15) 


p ^ 1 hnjkj 

’ MlkuT 


■]exp 


it) 


leq 


(16) 


where 


M = 


h{ajH + COC + 27r/) 

iEr ' 
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and are the empirical values associated with the observation of the 

magnetisation in the i,j G {x,\j,z\ directions relative to the hydrogen, carbon and 
J-coupled signals respectively. The time evolution of these three kinds of signals 
is presented schematically in Fig. In the figure one can see that at f = 0 , the 
hydrogen and carbon have the same magnetisation, as a consequence of the state 
preparation; at f = T^, the coherence has almost vanished; at f = T^, the system has 
a magnetisation reminiscent of the equilibrium, with the total hydrogen signal greater 
than that of the carbon; and at t = 60 s, which we call the system is practically 
thermalized, with the total signal equal to one, and the hydrogen signal almost 4 
times larger than that of the carbon. 



total signals intensities 


Figure 1. Comparison of the three distinct total signals intensities along time in the 
experiment. (Blue) Total hydrogen signals: = Y^i=x,y,z \ \ Eq- P^); (Yellow) 

total carbon signals: = Y^ i=x,y ,z I if (01 Eq. (15 T (Green) total coupled signals: 

I if (01 (cf- Eq. R). 


The equilibrium (Eq. ([13])) and non equilibrium states dealt with in the NMR 
experiment are of the kind 




1 + eA, 


(17) 


where e ~ 10 ^ and A is a traceless matrix, known as the deviation matrix, 
corresponding to the last three terms of Eq. (13). Though the experimental 


magnetisation signals come just from the deviation matrix, due to the traceless 
nature of the spin operators (for example: tr (pc; (8) Cj) = etr (Ac, (8) Cj)), these are 
macroscopic signals, for there are still a huge number ( 0 {e x 1023 ) = 0(101^)) of 
molecules contributing to the magnetisation. 

As discussed in the Introduction, the time evolution of the state simulated by 
a pseudo-pure state (Eq. (Q) is non-linear, if the process is non-unital. It is easy 
to circumvent this non-linearity. We have simply to tomograph the process acting 
directly on the NMR state given by Eq. ( [1^ , which of course evolves linearly. 

In summary, our strategy to keep measurements taken in different times properly 
normalized and comparable consists of two premisses, namely: (i) the quantum 
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tomography accounts for the process acting on the pseudo-pure state along the time, 
rather than the process acting on the simulated state; (ii) the process at any finite 
time must be compatible with the magnetisation converging to unity in an infinite 
time, i.e., the equilibrium magnetisation is taken as a reference to normalize all other 
measurements. 


4. The quantum tomographies 


The standard quantum process tomography (SQPT) [|^ consists in obtaining 
information about the map by means of quantum tomography of states modified by 
the action of the unknown process map. If the density matrices of the tomographed 
states have dimension d, then we need d^ linear independent states to reconstruct the 
quantum map. As we are dealing with two qubits, d is 4. Our experiment consists 
in preparing 20 different time average pseudo-pure states |j^| (q^, I = 1,2,... ,20), 
which can be used to simulate 20 different two-qubit pure states, where the density 
matrices of 16 of them are linear independent, spanning the Hilbert-Schmidt space. 
We also have 51 observations of these states in the time instants tj^ {k = 0 , 1 ,... , 50 ). 
The simulated 20 pure states correspond to 16 states forming the two-qubit Mutually 
Unbiased Basis (MUB) described in [20 1 , plus the four canonical states | tt)/ I ti)/ 

I I u) • 

The overcomplete set of 20 density matrices were preferred in order to increase 
the quality of the process tomography, following the arguments in [j^ [^. This 
strategy of enlarging the number of measurements with an overcomplete basis, in 
order to increase precision, was also used by one of us in the process tomography of 
a two-qubit map in Quantum Optics I30I. 

Since we have a large ensemble and sharp peaked measured signals 
we can consider the latter gaussianly distributed with standard deviation = 1 . 
Thus, we can deal with the experimental errors with a simple ^2"riorm minimization, 
which is equivalent to a log-likelihood method assuming gaussian noise 

The first thing to do is to perform a quantum state tomography in order to have 
a preliminary estimate of the initial states. We use a variation of ^ also equipped 
with nuclear norm minimization |}ii| as follows: 


min 

PO' 

subject to 


Po 

* k 

c,J} 



Po 

^ 0 




tr 


-lj'^(O) 

< 

2 

a''" 

tr 

^p[,l0C7f) 

1 

h 

< 

2 


tr 



< 

2 


Vi 

,j G {x,y,z} 





(18) 
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where || • ||^ stands for the nuclear norm or trace norm. This type of Semidefinite Problem 
(SDP) in the second-order cone can be efficiently solved using well known algorithms 

The program in Eq. ( [i^ has the following interpretation. We want a point 

estimate (pg) of a density matrix (positive semidefinite constraint), with minimum 

I i H C J\ 

residuals (A/y - ' in the measured data, according to minimization. The 

normalization of the states {pg} is unknown, and that is why it is convenient to adopt 
the trace norm minimization. The traces of the states {pg} are constrained by the 
normalization we imposed on the data, cf. Eq. ( [i3] > and Eq. ( [i^ , which corresponds 
to a unity magnetisation in equilibrium, cf. Eq. ( [i^ . 

After this preliminary estimation, we define a lower bound (Nq) for the trace of 
the states, namely: 


No = max < tr 


Po 


I = 1 ,..., 20 . 


(19) 


Then we rerun the program in Eq. (181 with the additional constraint: tr (pjj) > Nq. 


Note that the trace is related to the population in each state. The largest trace 
yielded by the first run of the state tomography identifies the population of responsive 
molecules. Therefore Nq is our reference population, and we are interested, at the initial 
time, in states with the same populations, in order to compare, along the time, the 
different expectation values. 

Now we consider the action of the unknown map ^{tk) upon the initial states 


determined in the second run of Eq. (181, namely: 




pk = ^{h){po} = E Xij{tk)APo^l 


(20) 


bi=i 


where the [Am] is a complete set of operators forming a basis in the Hilbert-Schmidt 
space. We obtain the basis {Am} by forming the linearly independent projectors 


corresponding to the states in our two-qubit MUB |2o|. The {Xijitk)} the elements 


of a positive semidefinite matrix, and form a representation of the quantum map in 
the basis {Am} @- 

We form the Liouville representation of the quantum state by putting the matrix 
elements of its d-dimensional density matrix, in lexicographical order, into a d^- 
dimensional vector jpg). We say jpg) is a "reshape" of pg. The Liouville representation 
of the map ^{t]^), in the operator basis {Am}, is a d^-dimensional matrix, namely Q: 

‘hfc = Y^Xijitk)Ai ® A*j. (21) 


Now the action of the map over the state is given by a simple matrix multiplication, 
namely: 


'Pk) =^k\po)- 


(22) 
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The density matrix |0^ is recovered by the reverse reshape operation. Now consider 
( 1 ^)} and {\tn)} as two orthonormal basis for a d-dimensional Hilbert space. If the 
matrix elements of are given by: 

{^k)nv = {mp\^j,\nv), (23) 

then the Choi matrix (or dynamical matrix) of the map is given by a "reshuffle" of 
<i>k, i.e., 

(Dt)»” = (24) 

For a map to be complete positive, its dynamical matrix must be positive semidefinite 
{Dk ^ 0)- A map is trace preserving if and only if 

UDk)’;," = h, ( 25 ) 

niji 

and it is unital if and only if 

E(D0™ = h- (26) 

nv 

At each time the quantum process tomography is obtained using a variation 
of the method in 1191, also equipped with a matrix recovery technique like nuclear 
norm minimization, in order to get the information about the trace: 

min ||Dt|L + EA''/,A''’ 

subject to xih) ^ 0 


Dfc A 0 


tr 


-y"(A) 

< A^'« 

2 ' 


tr 


u 

1 

< a'/C 

2 ^ 

(27) 

tr 


- il/M 

< 

2 - U 



Vi,y G {x,y,z} 

yi = i,..., 2 o 


The program in Eq. (27) can be interpreted as follows. We want a minimum trace 
positive semidefinite estimate of the dynamical matrix (Df), compatible with the data 
and having minimal residuals, analogously to the program of Eq. ( [1^ . In principle, 
it is redundant to include the constraint of positive semidefiniteness for both 
and xih)/ but, in practice, it gave us slightly better numerical precision using the 
interior-point based solvers MOSEK and SeDuMi [28 27I. 


5. Results and Discussion 

After running the program in Eq. ( [2^ , we plot in Eig. the trace distance between 
the states Pq imder the action of the map ^(t^) and the rescaled equilibrium state 



















H 



Figure 2. The relaxation of the states toward equilibrium. The figure shows the trace 
distance between the states pg under the action of the tomographed quantum map 
A(tj.) and the equilibrium thermal state. Note that the states in t = 0 are perturbations 
of the thermal state by different sequences of radiofrequency pulses. These perturbed 
thermal states are the ones useful for quanfum information processing simulafions. 


;^Peq- We note that the functional form of the curves are like an exponential decay, 
as expected. 

As a byproduct of the process tomography, the map must give us information 

CHI CHI 

about the relaxation times ' ' , T’2 ' ' . We define the quantity 





(28) 


where #l stands for the number of summed signals. This quantity is the average of 
the magnetisation in the f,/-direction {i,j G {x,y,z}) related to the equilibrium value 
and it gives us the functional behaviour of the magnetisation during the time. As a 
first attempt to model the functional behaviour of Eq. (28), one could fit a typical 
exponential decay like 


M(f) = + c, (29) 

where will be the aimed ^ This is of course a Markovian 

approximation of the xmderlying process. Notwithstanding, in the top panels of 
Figs, one can see slight deviations of the exponential decay. On the one hand, 
this could be a signature of non-Markovian behaviour, on the other hand, it could 
also be some artifact caused by the different state preparations necessary for each 
observation in time. We should emphasise that the tomography itself is independent 
of this exponential fit, but this simple model yields sensible results, in comparison 
with previously measured ' , T2' . The newly inferred T/, is also corroborated 
by Fig. which is independent of the exponential fitting. 

The results presented in Figs, [qjj^ and summarised in Table with the estimated 
errors, were computed using nonlinear least squares fitting with a trust-region 
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Figure 3. The Hydrogen {uncopled-)relaxation. (Top): the quantity (Mzi) as a function 

of time. (Bottom): the quantity function of time. Note that these 

"magnetisations" are relative to the equilibrium magnetisation, according to Eq. (281. 


algorithm pj] |. The least squares estimator was chosen to be consistent with the 
assumption that the experimental errors consist of uniform Gaussian noise, and also 
consistent with the £2“riorm minimization on the observed data in Eq. (271. 

The pseudo-pure states corresponding to the canonical states ( | tt)/ I ti)/ I it) / 
I it)) have magnetisation components (M21, Miz, Mzz) only in the z-direction. From 
these signals, we can retrieve the , as we can see in the top graphics in Figs. 
gand|^ 

The other 16 states have both longitudinal and transversal magnetisation 
components (M*i,Mi*, M**) and give us what we need to retrieve the T2' ' , as one 
can see in the bottom graphics in Figs. ^ and 

From the obtained results, one can conclude that, if one desires to simulate 
quantum parameters using the chloroform molecule, the time window where it can 
be done is 0 . 177 s, which is the lifetime of the coherent coupling of the spins (T^) in 
the chloroform molecule. Otherwise, we can consider the system (partially) relaxed 
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Figure 4. The Carbon {uncoupled-)relaxation. (Top): the quantity (Mjz) as a function 
of time. (Bottom): the quantity ^ function of time. Note that these 


"magnetisations" are relative to the equilibrium magnetisation, according to Eq. (281. 



ri(s) 

Tlis) 

Hydrogen 

6.2 ± 0.2 

0.24 ± 0.01 

Carbon 

7.4 ± 0.3 

0.192 ± 0.005 

J-coupling 

5.65 ± 0.07 

0.177 ± 0.001 


Table 1. Relaxation times retrieved from the curve fit in the Figs. and|^using the 
model in Eq. |[^. 


and unable to simulate all parameters properly for a quantum computation. In 
Fig ^ we confirm that is the characteristic time window for quantum computation 
simulations by means of the properties of the quantum process relaxation map. One 
clearly see in that figure that the map is almost perfectly unital and trace-preserving 
up to log(r2 = 0 . 177 s) = — 1 . 731 . After rj the map dramatically departs from been 
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Figure 5. The J-coupling relaxation. (Top): the quantity (Mzz) as a funchon of time. 

(Bottom): the quantity + ^ function 

of time. Notice how small is this coupled magnetisation compared to the uncoupled 
magnetisations shown in Figs. 3 and 4. Note also that these "magnetisations" are 
relative to the equilibrium magnetisation, according to Eq. |[2^. 


both trace-preserving and unital, and therefore the normalization problem comes 
into play. These results demonstrate the robustness and flexibility of our quantum 
tomography methodology 

6. Conclusion 

Summarising, we performed a quantum process tomography of the relaxation of 
a two-spin system in an NMR liquid-state experiment. As the experimental data 
correspond to expectation values of traceless operators, the use of informational 
incomplete quantum tomography techniques to handle the missing quantum state 
normalization information was paramoimt to reconstruct the time dependent non¬ 
uni tal and non trace-preserving quantum map. The successful approach we employed 
allowed us to introduce two new parameters (T|, T^) that are relevant for the 
characterization of the relaxation process. We learned that the coupled magnetisations 
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log( time (s)) 


Figure 6. Quantification of how much the relaxation process departs from being 
unital {cf. Eq. | |26^ ) (black circles), and trace-preserving {cf. Eq. | |25^ ) (red triangles). 
The dramatic change in the behaviour of the map occurs after log ~ — 1 . 8 . 

have a similar behaviour to the uncoupled ones, but with different characteristic 
times. The quantum process tomography showed that characterizes the time 
window for the simulation of quantum computation. It can be seen clearly in Figj^ 
For the simulation of quantum maps for times greater than T^, one has to accoimt 
for the time dependent normalization of the data, for the map associated to the 
ever present relaxation process cannot be assumed to be unital and trace-preserving. 
We also introduced an heuristic for the normalization problem, which consists in 
assuming the theoretical equilibrium magnetisation as the reference to compare 
all the measurements, and the largest trace obtained with the preliminary state 
tomography (Eq. ( [i8] )) as a lower bound for the trace of all states. 

7. Acknowledgements 

We acknowledge Johan Aberg, Paulo H. S. Ribeiro, Lucas C. Celeri and Jefferson G. 
Filgueiras for fruitful discussions. We are indebted to the anonymous referees for their 
criticism and suggestions that largely contributed to the final version of this paper. 
Financial support by Brazilian agencies CAPES, CNPq, FAPERJ, EAPEMIG, and 
INCT-IQ (National Institute of Science and Technology for Quantum Information). 







REFERENCES 


19 


References 

[1] Koenraad M R Audenaert and Stefan Scheel. 'Quantum tomographic 
reconstruction with error bars: a Kalman filter approach'. In: New Journal 
of Physics 11.2 (2009), p. 023028. doi: 10 . 1088 / 1367 - 2630 / 11 / 2 / 023028 , 

[2] 1 . Bengtsson and K. Zyczkowski. Geometry of Quantum States, An Introduction to 
QUANTUM ENTANGEEMENT. 1st ed. Cambridge University Press, 2008. 

[3] N Boulant, T F Havel, M A Pravia and D G Cory. 'Robust method for estimating 
the Lindblad operators of a dissipative quantum process from measurements 
of the density operator at multiple time points'. In: Physical Review A 67.4 (Apr. 
2003), pp. 042322-12. 

[4] S. Boyd and L. Vandenberghe. Convex Optimization. Cambridge University Press, 
2004. 

[5] Mark D. de Burgh, Nathan K. Langford, Andrew C. Doherty and Alexei 
Gilchrist. 'Choice of measurement sets in qubit tomography'. In: Phys. Rev. 
A 78 (5 2008), p. 052122. doi: 10 . 1103 /PhysRevA . 78 . 052122 . url: http : 
//link.aps.org/doi/ 10 . 1103 /PhysRevA. 78 . 052122 , 

[6] Isaac L. Chuang and M. A. Nielsen. 'Prescription for experimental 
determination of the dynamics of a quantum black box'. In: Journal of Modern 
Optics 44.11-12 (1997), pp. 2455-2467. doi: 10.1080/09500349708231894, 

[7] D G Cory, A F Fahmy and T F Havel. 'Ensemble quantum computing by NMR 
spectroscopy'. In: Proceedings of the National Academy of Sciences of the United 
States of America 94.5 (1997), pp. 1634-1639. 

[8] A Gavini-Viana, A M Souza, D O Soares-Pinto, J Teles, R S Sarthour, E R 
deAzevedo, T J Bonagamba and 1 S Oliveira. 'Normalization procedure for 
relaxation studies in NMR quantum information processing'. In: Quantum 
Information Processing 9.5 (Dec. 2009), pp. 575-589. 

[9] D. S. Gonsalves, C. Lavor, M. A. Gomes-Ruggiero, A. T. Cesario, R. O. Vianna 
and T. O. Maciel. 'Quantum state tomography with incomplete data: Maximum 
entropy and variational quantum tomography'. In: Phys. Rev. A 87 (5 2013), 
p. 052140. doi: 10 . 1103 /PhysRevA. 87 . 052140 . 

[10] Christopher E Granade, Christopher Eerrie, Nathan Wiebe and D G Cory. 
'Robust online Hamiltonian learning'. In: New Journal of Physics 14.10 (Oct. 
2012), pp. 103013-32. 

[11] David Gross, Yi-Kai Liu, Steven T. Elammia, Stephen Becker and Jens Eisert. 
'Quantum State Tomography via Compressed Sensing'. In: Phys. Rev. Eett. 105 
(15 2010), p. 150401. doi: 10 . 1103 /PhysRevLett. 105 . 150401 , 

[12] E. T. Jaynes. 'Information Theory and Statistical Mechanics'. In: Physical Review 
106.4 (1957)/ PP- 620-630. 



REFERENCES 


20 


[13] T. D. Ladd, F. Jelezko, R. Laflamme, Y. Nakamura, C. Monroe and J. L. O'Brien. 
'Quantum computers'. In: Nature 464.7285 (2010), pp. 45-53. url: http://dx 
doi.org/ 10 . 1038 /nature 08812 , 

[14] Nathan K Langford. 'Errors in quantum tomography: diagnosing systematic 
versus statistical errors'. In: New Journal of Physics 15.3 (2013), p. 035003. url: 
http://stacks.iop.org/ 1367 - 2630 / 15 / 1 = 3 /a= 035003 , 

[15] Garett M. Leskowitz and Leonard J. Mueller. 'State interrogation in nuclear 
magnetic resonance quantum-information processing'. In: Phys. Rev. A 69 (5 
2004), p. 052302. doi: 10 . 1103 /PhysRevA. 69 . 052302 , 

[16] G. Lima, E. S. Gomez, A. Vargas, R. O. Vianna and C. Saavedra. 'East 
entanglement detection for unknown states of two spatial qutrits'. In: Phys. 
Rev. A 82 (1 2010), p. 012302. doi: 10 . 1103 /PhysRevA. 82 . 012302 , url: http: 
//link.aps.org/doi/ 10 . 1103 /PhysRevA. 82 . 012302 , 

[17] J. Lofberg. YAEMIP : A Toolbox for Modeling and Optimization in MATEAB. Taipei, 
Taiwan, 2004. url: http://control. ee.ethz.ch/\~{}-joloef/yalmip.php, 

[18] T. O. Maciel, A. T. Cesario and R. O. Vianna. 'Variational quantum 
tomography with incomplete information by means of semidefinite programs'. 
In: International Journal of Modern Physics C 22.12 (2011), pp. 1361-1372. 

[19] T. O. Maciel and R. O. Vianna. 'Optimal estimation of quantum processes 
using incomplete information: variational quantum process tomography'. In: 
Quantum Information and Computation 12.5 (2012), pp. 0442-0447. 

[20] T. O. Maciel and R. O. Vianna. 'Viable entanglement detection of unknown 
mixed states in low dimensions'. In: Phys. Rev. A 80 (3 2009), p. 032325. doi: 
10 . 1103 /PhysRevA. 80 . 032325 , 

[21] MATLAB. version 8.y.o.igy6i^ (Rzoiya). Natick, Massachusetts: The Math Works 
Inc., 2015. 

[22] M. Nielsen and 1 . Chuang. Quantum Computation and Quantum Information. 
1st ed. Cambridge University Press, 2000. 

[23] M. A. Nielsen, E. Knill and R. Laflamme. 'Complete quantum teleportation 
using nuclear magnetic resonance'. In: Nature 396.6706 (1998), pp. 52-55. url: 
http://dx.doi.org/ 10 . 1038 / 23891 . 

[24] 1 . Oliveira, R. Sarthour, T. Bonagamba, E. Azevedo and J.C.C. Ereitas. NMR 
Quantum Information Processing. Elsevier Science, 2011. 

[25] W M Pimenta, B Marques, T O Maciel, R O Vianna, A Delgado, C Saavedra 
and S Padua. 'Minimum tomography of two entangled qutrits using local 
measurements of one-qutrit symmetric informationally complete positive 
operator-valued measure'. In: Physical Review A 88.1 (July 2013), pp. 012112-8. 

[26] M Ringbauer, C J Wood, K Modi, A Gilchrist, A G White and A 
Pedrizzi. 'Characterizing Quantum Dynamics with Initial System-Environment 
Correlations'. In: Physical Review Fetters 114.9 (Mar. 2015), pp. 090402-5. 



REFERENCES 


21 


[27] Jos R Sturm. Using SeDuMi 1.02, a MATEAB toolbox for optimization over 
symmetric cones. 1998. 

[28] The MOSEK optimization software, http://www.mosek.com/. url: http : //www 
mosek.com/. 

[29] Lieven M. K. Vandersypen, Matthias Steffen, Gregory Breyta, Costantmo S. 
Yannoni, Mark H. Sherwood and Isaac L. Chuang. 'Experimental realization 
of Shor's quantum factoring algorithm using nuclear magnetic resonance'. In: 
Nature 414.6866 (Dec. 2001), pp. 883-887. url: http : //dx. doi . org/ 10 . 1038 / 
414883 a, 

[30] Reinaldo O Vianna, Andrea Crespi, Roberta Ramponi, Roberto Osellame, Linda 
Sansoni, Giorgio Milani, Paolo Mataloni and Fabio Sciarrino. 'Variational 
quantum process tomography of two-qubit maps'. In: Physical Review A 87.3 
(2013). 

[31] Yaakov S Weinstein, Timothy F Havel, Joseph Emerson, Nicolas Boulant, Marcos 
Saraceno, Seth Lloyd and David G Cory. 'Quantum process tomography of the 
quantum Fourier transform'. In: The Journal of Chemical Physics 121.13 (2004), 
pp. 6117-18. 




